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Abstract 

Using Dyson-Schwinger equations we obtain an ultraviolet asymptotics for the dynamical 
mass of quark in QCD 



We also determine a numerical value for the vr meson decay constant fj^. 
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As is well known, strong interactions at high energies are approximately invariant with respect 
to the chiral symmetry group SUl{N) (g) SUr{N) {N is the number of different quark species). 
The Lagrangian of QCD - the most popular field-theoretical model of strong interactions - has 
the chiral symmetry group SUl{N) (g) SUfi,{N) in approximation of massless quarks. The success 
of current algebra and PCAC hypothesis shows that the chiral symmetry SUl{N) ® SUr{N) 
is spontaneously broken down to SU{N) symmetry. The spontaneous violation of the chiral 
symmetry reveals itself in nonvanishing vacuum expectation values (qq) and also in the appearance 
of dynamical mass for originally massless quarks. 

The problem of dynamical mass generation in massless theories has a long history. First, 
in refs. [0, |^ it has been shown in some concrete models that the generation of the dynamical 
mass can take place in quantum field theory. However, the consistent study of the problem of 
dynamical mass generation requires going beyond perturbation theory. It is obvious that this 
problem is rather nontrivial in the field theory because of the absence of the consistent methods 
different from perturbation theory. 

The purpose of the present paper is to investigate a possibility of generation of the dynamical 
mass of quarks in QCD. Using the Dyson- Schwinger equations we show that one can reliably 
determine the dependence of the dynamical mass of quarks m{p) on the momentum p in the 
ultraviolet domain due to the property of asymptotic freedom of strong interactions. In the 
following we will work in Euclidean space-time and in the transverse gauge (Landau gauge). The 
use of this gauge is convenient since the propagator of the massless quark is not renormalized in 
the leading logarithmic approximation. 

The Dyson-Schwinger equation for the quark propagator G{p) has the form 
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G-\p) = S-\p) + ^/ d'kT^^{p,k)G{k)D,,{p- k)^,{Xy2). (1) 

Here S{p) is the propagator of free quarks, A"/2 are the generators of the fundamental represen- 
tation of the group SU{N). Equation (0) is written for bare quantities. After performing the 
renormalization we find 
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Z2'Gj,\p) = S-\p) + ^/ d'kTf{p,k)Gn{k)D^,{p-k)^^{\y2), 

where Z2 is the renormalization constant of the quark propagator. At the lowest order in gj^ in the 
Landau gauge the renormalization constant of the quark propagator is equal to one Z2 = 1+0 {gj^). 
Therefore in the leading logarithmic approximation we can set G^^{p) = m{p)—p. For the quantity 
m{p) we obtain the equation of the form 



m{p ) = d k 



m{k'^)g'^{{p — k) 



2^ 



{2tt)^J {m?{k^) + k^){p-kY 

where 

f{p') = 9\v'l^\9l) = 91/ (1 + (11 - ^m^2^^^2 

is the effective quark-gluon coupling constant. We split the region of integration into two sub- 
domains determined by the requirements k'^ > p^ (sub-domain (I)) and k"^ < p^ (sub-domain (II)) 

d^k= I d^k+ I d^k. 

Jk^yp-^, (/) Jk^<p2, {II) 
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Then, taking into account that 



d _2, s 33 2N 2/ \\2 
(-) = -^^(^ (-)) 



we obtain the equahty g'^[{p — kY\ = g'^{k'^) up to the terms of the next order in in the sub- 
domain (I) where /c^ > and the equahty g^[{p — /c)^] = g'^ip^) with the same accuracy in the 
sub-domain (II) where > k"^. 

Setting = X we obtain after integrating over the angular variables 

1 g'^{z)m{z)dz 1 g^{x) Z"^' m{z)zdz , , 

"^W = 4;^y^ ,„2(^) + ^ +47,2 ^ m^{z) + z' 

After some transformations of eq. (0) one finds 

d I Att"^ {dm / dx)\ m{x)x 



dx \ g'^ix) J 171"^ (x) + X 

Taking into account that 

(33 - 2N)f{x) = A8ny\n{x/A^) 
and changing the variable x/A^ = exp(t) we obtain 

d'^T.(t) , dY.it) Yit) , 9.^. 9 . 1 

+ (1 + lA)^ + + S2(t)A-2 exp(-t))-i = (3) 

where d = 12/(33 - 2N), m{x) = S(t) with t = ln(x/A2). 

We are looking for a solution of eq. (^) which is bounded at t — > cxo, therefore in the asymptotic 
regime one can neglect the nonlinearity of the equation. The linear equation has the solution 

S(t) = exp(-t/2 - ^ Int) (C, W.^+iM't) + C2 W^d--i/2,o(t) 
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where W\_fj_{z) are Whittaker functions with the asymptotic behavior for large z 

W^A,^(^)||.|-.oo = exp(-^/2)z^[l + 0(1/^)] . (4) 
Using eq. (^ we obtain the asymptotics of two linearly independent solutions of eq. (^: 

mi(x)|™~Ciln-'^(a;/A2) 

and 

m2(x)|™~(C2/x) ln'^-^(x/A2). 

For the determination of the constants Ci and C2 we use the comparison with perturbation 
theory. We require smallness of the nonperturbative corrections in the asymptotic domain in the 
form 

m{p^) — fh{p^) = d{fh{p^)), ^ 00 

where fh(p'^) = Tn^ijj'^ijP') /g'^^ji^))''' is the usual expression for the effective quark mass in the 
renormalization group analysis. Then 
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and 

m{p') = mn {fip') I f {^^'))' - j,f{p') [fip') I f {l^'))'" ■ (5) 

The second term in this expression describes the ultraviolet asymptotics of the dynamical mass 
of the light quark. For the determination of the constant C we use the perturbation theory with 
phenomenological account for the nonperturbative effects through (gg) 7^ 0. 

In the second order of the perturbation theory with the nonvanishing vacuum expectation 
value of the bilinear quark operator {qq) one has 

on the other hand the use of the explicit form of the mass operator from eq. gives 

tr G{p) = -ACgl/{p'f. 

Comparing these two expressions we obtain 

C = (gg)/3. (6) 

With the expression for the dynamical mass given in eq. where the constant C is fixed by 
the relation from eq. (BI) we compute the pion decay constant /^r defined through 



where |7r"(/c)) is the one-particle pion state, A^(0) is the isotopic vector axial current. The rep- 
resentation for the pion decay constant /tt was found in ref. P] and used earlier for a similar 
calculation in ref. @] 

3 roo X {m{x) - ^'^) m{x) 



(27r)2 JO {x + m'^{x)y 

For determining the numerical values one has to know the magnitude of the vacuum expectation 
value (qq). The methods of current algebra allow one to relate (qq) to the current masses of light 
quarks that leads to the estimate (qq) = —(250 MeV)^. Then the relation from eq. (0) gives 
/tt = 77 MeV. Using the experimental value /tt = 93 MeV one can turn the problem around 
and determine \ {qq) \ = (310 MeV)^. Note that in this way one cannot determine the sign of the 
vacuum expectation value (qq). 
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